Next-generation RNA Sequencing Research Articles

Abstract PO1-14-08: Real-world clinical genomics study of HR+/HER2- metastatic breast cancers treated by CDK4/6i plus endocrine therapies revealed a drug resistant tumor segment characterized by ER independence

Abstract Background CDK4/6 inhibitors (CDK4/6i) plus endocrine therapies (ET) are the standard-of-care for hormone receptor–positive/human epidermal receptor 2–negative metastatic breast cancer (HR+/HER2− mBC). However, drug resistance remains a major unmet need. Investigations of drug resistance mechanisms has been hampered by a dearth of tumor molecular profiling data from the post-treatment setting. To address this challenge, we have conducted a real-world clinical genomics study to better understand the molecular mechanism of CDK4/6i resistance as well as to stratify patients based on integrated multi-omics profiles. Methods We retrospectively analyzed a multi-omics dataset of 400 HR+/HER2- mBC patients who had received CDK4/6i plus ET and developed progressive disease (PD) from the de-identified Tempus database. Pre-treatment and post-progression biopsies were taken  1 year prior to starting the CDK4/6i treatment or following PD respectively. Tempus xT next-generation sequencing (DNA-seq of 648 genes) and RNA sequencing assays were performed on 427 tumor FFPE samples, including 200 pre-treatment, 227 post-progression and 26 longitudinal pairs. Results The median age of the patients was 57 (54.9-57.4) and median progression free survival (PFS) is 379 (341-433) days. Two genes were found to harbor a significant increase in genomic alteration frequencies (GAF) after adjusting for FDR at post-progression vs. pre-treatment – ESR1 (41.9% vs. 15%, p=5.4e-10), RB1 (13.2% vs. 3%, p=8.5e-05). ESR1 and RB1 also harbored high frequencies of acquired genomic alterations among 26 paired samples at 34.6% and 11.5% respectively. TP53 mutation at baseline was significantly associated with shorter PFS at baseline (p=4.23e-05, HR=2.081) and TP53 GAF significantly increased after PD (37% vs. 28.5%, p=0.039). BRCA1/2 pathogenic mutations (p=1.63e-04, HR=3.066), APOBEC mutation signature S13 (p=0.0125, HR=1.55) and CCNE1 gene expression (p=0.024, HR=1.46) were significantly associated with shorter PFS. APOBEC signature (p=0.0035) and CCNE1 expression (p=1.33e-06) also significantly increased post-progression. Among the top molecular features associated with longer PFS were markers of estrogen signaling such as PGR gene expression (p=6.76e-04, HR=0.565) and the Hallmark estrogen response signature (p=0.021, HR=0.679). Applying a multi-omics pattern recognition algorithm, we identified a molecularly distinct cluster (IC1) characterized by down-regulation of estrogen signaling. IC1 is significantly associated with shorter PFS (p=3.72e-05, HR=0.22) and increased from 4% pre-treatment to 23% post-progression (p=7.3e-08). Further, IC1 is strongly enriched in markers previously implicated in CDK4/6i resistance including CCNE1 expression, RB1 mutation and MYC/E2F activation. We then developed machine learning models to predict gene-level dependency trained on cancer cell line expression and CRISPR-KO screen data. These models predicted decreased dependency on ESR1 and CDK4 and increased dependency on CDK2 in IC1, strengthening the association between ER independence and CDK4/6i resistance. Conclusions Our real-world clinical genomics study identified a comprehensive list of biomarkers associated with resistance to CDK4/6i plus ET and estimated patient prevalence for these markers in the post-treatment setting. Integrated and machine-learning analyses identified a subset of aggressive tumors with estrogen independence characteristics that are implicated in CDK4/6i resistance and suggested new therapeutic strategies. Citation Format: Zhengyan Kan, Ji Wen, Jennifer Webster, Vinicius Bonato, Whijae Roh, Xinmeng Jasmine Mu, Paul Rejto, Jadwiga Bienkowska. Real-world clinical genomics study of HR+/HER2- metastatic breast cancers treated by CDK4/6i plus endocrine therapies revealed a drug resistant tumor segment characterized by ER independence [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PO1-14-08.

Genomic Profiling of Rare Undifferentiated Sarcomatoid Subtypes of Pancreatic Carcinomas: In Search of Therapeutic Targets

PURPOSE The highly aggressive undifferentiated sarcomatoid carcinoma (USC) subtype of pancreatic ductal adenocarcinoma (PDAC) remains poorly characterized because of its rarity. Previous case reports suggest that immune checkpoint inhibitors could be a promising treatment strategy, but the prevalence of established predictive biomarkers of response is largely unknown. The objective of this study was to leverage comprehensive genomic profiling of USC PDAC tumors to determine the prevalence of biomarkers associated with potential response to targeted therapies. METHODS USC tumors (n = 20) underwent central pathology review by a board-certified gastrointestinal pathologist to confirm the diagnosis. These samples were compared with non-USC PDAC tumors (N = 5,562). Retrospective analysis of DNA and RNA next-generation sequencing data was performed. RESULTS USC PDACs were more frequently PD-L1+ by immunohistochemistry than non-USC PDAC (63% v 16%, respectively, P < .001). Furthermore, USC PDAC had an increase in neutrophils (8.99% v 5.55%, P = .005) and dendritic cells (1.08% v 0.00%, q = 0.022) and an increased expression of PDCD1LG2 (4.6% v 1.3%, q = 0.001), PDCD1 (2.0% v 0.8%, q = 0.060), and HAVCR2 (45.9% v 21.7%, q = 0.107) than non-USC PDAC. Similar to non-USC PDAC, KRAS was the most commonly mutated gene (86% v 90%, respectively, P = 1). CONCLUSION To our knowledge, this work represents the largest molecular analysis of USC tumors to date and showed an increased expression of immune checkpoint genes in USC tumors. These findings provide evidence for further investigation into immune checkpoint inhibitors in USC tumors.

Cytosolic DNA initiates a vicious circle of aging-related endothelial inflammation and mitochondrial dysfunction via STING: the inhibitory effect of Cilostazol.

Endothelial senescence, aging-related inflammation, and mitochondrial dysfunction are prominent features of vascular aging and contribute to the development of aging-associated vascular disease. Accumulating evidence indicates that DNA damage occurs in aging vascular cells, especially in endothelial cells (ECs). However, the mechanism of EC senescence has not been completely elucidated, and so far, there is no specific drug in the clinic to treat EC senescence and vascular aging. Here we show that various aging stimuli induce nuclear DNA and mitochondrial damage in ECs, thus facilitating the release of cytoplasmic free DNA (cfDNA), which activates the DNA-sensing adapter protein STING. STING activation led to a senescence-associated secretory phenotype (SASP), thereby releasing pro-aging cytokines and cfDNA to further exacerbate mitochondrial damage and EC senescence, thus forming a vicious circle, all of which can be suppressed by STING knockdown or inhibition. Using next-generation RNA sequencing, we demonstrate that STING activation stimulates, whereas STING inhibition disrupts pathways associated with cell senescence and SASP. In vivo studies unravel that endothelial-specific Sting deficiency alleviates aging-related endothelial inflammation and mitochondrial dysfunction and prevents the development of atherosclerosis in mice. By screening FDA-approved vasoprotective drugs, we identified Cilostazol as a new STING inhibitor that attenuates aging-related endothelial inflammation both in vitro and in vivo. We demonstrated that Cilostazol significantly inhibited STING translocation from the ER to the Golgi apparatus during STING activation by targeting S162 and S243 residues of STING. These results disclose the deleterious effects of a cfDNA-STING-SASP-cfDNA vicious circle on EC senescence and atherogenesis and suggest that the STING pathway is a promising therapeutic target for vascular aging-related diseases. A proposed model illustrates the central role of STING in mediating a vicious circle of cfDNA-STING-SASP-cfDNA to aggravate age-related endothelial inflammation and mitochondrial damage.

MiRNome profiling analysis reveals novel hepatocellular carcinoma diagnostic, prognostic and treatment-related candidate biomarkers: post-hoc analysis of SORAMIC trial.

Introduction Early diagnosis of hepatocellular carcinoma as well as evaluation of prognosis and prediction of treatment efficacy remain challenging due to the missing specific non-invasive biomarkers. The aim of this study is to identify disease-specific microRNA (miRNA) patterns for diagnosis, prediction of prognosis and treatment response in patients with hepatocellular carcinoma (HCC). Methods The study population included 42 HCC patients from SORAMIC clinical trial: 22 patients received sorafenib monotherapy, 20 patients underwent 90Y radioembolization in combination with sorafenib. 20 individuals were included in the control group. Hepatocellular carcinoma patients underwent collection of plasma samples before and 7-9 weeks after the beginning of the treatment. Isolation of circulating miRNAs, preparation of small RNA sequencing libraries and next-generation sequencing were performed. Association analysis for novel diagnostic, prognostic and treatment-related candidate biomarkers was performed. Results A total of 42 differentially expressed (16 up-regulated and 26 down-regulated) miRNAs were identified comparing baseline and control group plasma samples. hsa-miR-215-5p and hsa-miR-192-5p were down-regulated, while hsa-miR-483-5p and hsa-miR-23b-3p were up-regulated comparing baseline and 7-9 weeks post-sorafenib monotherapy samples. hsa-miR-215-5p was the sole down-regulated miRNA in the same combination therapy comparison. hsa-miR-183-5p, hsa-miR-28-3p and hsa-miR-1246 were found to be significantly up-regulated comparing non-responders versus responders to sorafenib. High hsa-miR-215-5p expression was significantly associated with worse HCC patients' prognosis. Conclusions Systematic miRNA profiling of highly characterized samples from SORAMIC study revealed a subset of potential miRNA biomarkers for hepatocellular carcinoma diagnosis and prognosis of sorafenib-treated patients' survival.

CD58 alterations govern antitumor immune responses by inducing PD-L1 and IDO in diffuse large B-cell lymphoma.

Recurrent abnormalities in immune surveillance-related genes affect the progression of diffuse large B-cell lymphoma (DLBCL) and modulate the response to therapeutic interventions. CD58 interacts with the CD2 receptor on T cells and natural killer (NK) cells and is recurrently mutated and deleted in DLBCL, suggesting it may play a role in regulating antitumor immunity. Herein, we comprehensively analyzed the genomic characteristics of CD58 through targeted next-generation sequencing, RNA-sequencing, whole-exome sequencing, and single-cell RNA-sequencing in patients with newly diagnosed DLBCL. The CD58 mutation rate was 9.1%, and the copy number loss rate was 44.7% among all enrolled DLBCL patients. Notably, CD58 genetic alterations, along with low CD58 expression, significantly correlated with reduced rates of response to R-CHOP therapy and inferior progression-free and overall survival. Single-cell RNA sequencing revealed that CD58 expression in tumor cells was negatively correlated with CD8+ T cell exhaustion/dysfunction status. Insufficient T-cell activation resulting from CD58 alterations could not be attributed solely to CD2 signaling. CD58 inhibited the activity of the JAK2/STAT1 pathway by activating the Lyn/CD22/SHP1 axis, thereby limiting PD-L1 and IDO expression. Elevated PD-L1 and IDO expression in CD58 deficient DLBCL cells led to immune evasion and tumor-intrinsic resistance to CAR T-cell therapy. Direct activation of CD58-CD2 costimulatory signaling in combination with anti-PD-L1 blockade or IDO inhibitor sensitized CD58-deficient DLBCL to CAR T-cell therapy. Collectively, this work identified the multiple roles of CD58 in regulating antitumor immune responses in DLBCL.

A comprehensive analysis of pneumococcal two-component system regulatory networks.

Two-component systems are key signal-transduction systems that enable bacteria to respond to a wide variety of environmental stimuli. The human pathogen, Streptococcus pneumoniae (pneumococcus) encodes 13 two-component systems and a single orphan response regulator, most of which are significant for pneumococcal pathogenicity. Mapping the regulatory networks governed by these systems is key to understand pneumococcal host adaptation. Here we employ a novel bioinformatic approach to predict the regulons of each two-component system based on publicly available whole-genome sequencing data. By employing pangenome-wide association studies (panGWAS) to predict genotype-genotype associations for each two-component system, we predicted regulon genes of 11 of the pneumococcal two-component systems. Through validation via next-generation RNA-sequencing on response regulator overexpression mutants, several top candidate genes predicted by the panGWAS analysis were confirmed as regulon genes. The present study presents novel details on multiple pneumococcal two-component systems, including an expansion of regulons, identification of candidate response regulator binding motifs, and identification of candidate response regulator-regulated small non-coding RNAs. We also demonstrate a use for panGWAS as a complementary tool in target gene identification via identification of genotype-to-genotype links. Expanding our knowledge on two-component systems in pathogens is crucial to understanding how these bacteria sense and respond to their host environment, which could prove useful in future drug development.

Insights from bioinformatics analysis reveal that lipopolysaccharide induces activation of chemokine-related signaling pathways in human nasal epithelial cells

Lipopolysaccharide (LPS) is known to elicit a robust immune response. This study aimed to investigate the impact of LPS on the transcriptome of human nasal epithelial cells (HNEpC). HNEpC were cultured and stimulated with LPS (1 μg/mL) or an equivalent amount of normal culture medium. Subsequently, total RNA was extracted, purified, and sequenced using next-generation RNA sequencing technology. Differentially expressed genes (DEGs) were identified and subjected to functional enrichment analysis. A protein–protein interaction (PPI) network of DEGs was constructed, followed by Ingenuity Pathway Analysis (IPA) to identify molecular pathways influenced by LPS exposure on HNEpC. Validation of key genes was performed using quantitative real-time PCR (qRT-PCR). A total of 97 DEGs, comprising 48 up-regulated genes and 49 down-regulated genes, were identified. Results from functional enrichment analysis, PPI, and IPA indicated that DEGs were predominantly enriched in chemokine-related signaling pathways. Subsequent qRT-PCR validation demonstrated significant upregulation of key genes in these pathways in LPS-treated HNEpC compared to control cells. In conclusion, LPS intervention profoundly altered the transcriptome of HNEpC, potentially exacerbating inflammatory responses through the activation of chemokine-related signaling pathways.

The plasma miRNome and venous thromboembolism in high-grade glioma: miRNA Sequencing of a nested case-control cohort.

Patients with high-grade gliomas are at high risk of venous thromboembolism (VTE). MicroRNAs (miRNAs) are small non-coding RNAs with multiple roles in tumour biology, haemostasis and platelet function. Their association with VTE risk in high-grade glioma has not been comprehensively mapped so far. We thus conducted a nested case-control study within 152 patients with WHO grade IV glioma that had been part of a prospective cohort study on VTE risk factors. At inclusion a single blood draw was taken, and patients were thereafter followed for a maximum of 2 years. During that time, 24 patients (16%) developed VTE. Of the other 128 patients, we randomly selected 24 age- and sex-matched controls. After quality control, the final group size was 21 patients with VTE during follow-up and 23 without VTE. Small RNA next-generation sequencing of plasma was performed. We observed that hsa-miR-451a was globally the most abundant miRNA. Notably, 51% of all miRNAs showed a correlation with platelet count. The analysis of miRNAs differentially regulated in VTE patients-with and without platelet adjustment-identified potential VTE biomarker candidates such as has-miR-221-3p. Therewith, we here provide one of the largest and deepest peripheral blood miRNA datasets of high-grade glioma patients so far, in which we identified first VTE biomarker candidates that can serve as the starting point for future research.

Oral and gut microbial profiling in periodontitis and Parkinson’s disease

ABSTRACT Objectives We tested the hypothesis that Parkinson’s disease (PA) alters the periodontitis-associated oral microbiome. Method Patients with periodontitis with Parkinson’s disease (PA+P) and without PA (P) and systemically and periodontally healthy individuals (HC) were enrolled. Clinical, periodontal and neurological parameters were recorded. The severity of PA motor functions was measured. Unstimulated saliva samples and stool samples were collected. Next-generation sequencing of 16S ribosomal RNA (V1-V3 regions) was performed. Results PA patients had mild-to-moderate motor dysfunction and comparable plaque scores as those without, indicating that oral hygiene was efficient in the PA+P group. In saliva, there were statistically significant differences in beta diversity between HC and PA+P (p = 0.001), HC and P (p = 0.001), and P and PA+P (p = 0.028). The microbial profiles of saliva and fecal samples were distinct. Mycoplasma faucium, Tannerella forsythia, Parvimonas micra, and Saccharibacteria (TM7) were increased in P; Prevotella pallens, Prevotella melaninogenica, Neisseria multispecies were more abundant in PA+P group, Ruthenibacterium lactatiformans, Dialister succinatiphilus, Butyrivibrio crossotus and Alloprevotella tannerae were detected in fecal samples in P groups compared to healthy controls. Conclusions No significant differences were detected between Parkinson’s and non-Parkinson’s gut microbiomes, suggesting that Parkinson’s disease modifies the oral microbiome in periodontitis subjects independent of the gut microbiome.

Next-generation sequencing-based detection reveals Erysiphe necator-associated virus 1 in okra plants

Erysiphe necator is a devastating fungal phytopathogen that has caused powdery mildew on several genera within the Vitaceae family worldwide. In the current investigation, next-generation sequencing of RNA extracted from the okra crop, besides the bioinformatics analysis, was operated to examine the virome associated with the okra crop. Although E. necator fungus was not identified in the okra tissues, the mycovirus-related sequence of Erysiphe necator-associated virus 1 was detected. This mycovirus sequence showed a high similarity ranging between 80.38-93.48% and 100% coverage with the four distinct global isolates of the Erysiphe necator-associated virus 1. This genomic data postulates comprehension of the existence of the mycoviruses in plant hosts without their primary host, the fungus. To our knowledge, this is the first identification of Erysiphe necator-associated virus 1 in Iraq

Bulk and Single-Cell RNA Sequencing Elucidate the Etiology of Severe COVID-19.

Coronavirus disease 2019 (COVID-19) is a type of pneumonia caused by the SARS-CoV-2 coronavirus. It can cause acute pulmonary and systemic inflammation, which can lead to death in severely ill patients. This study explores the potential reasons behind severe COVID-19 and its similarities to systemic autoimmune diseases. This study reviewed unbiased high-throughput gene expression datasets, including next-generation and single-cell RNA sequencing. A total of 27 studies and eight meta-analyses were reviewed. The studies indicated that severe COVID-19 is associated with the upregulation of genes involved in pro-inflammatory, interferon, and cytokine/chemokine pathways. Additionally, changes were observed in the proportions of immune cell types in the blood and tissues, along with degenerative alterations in lung epithelial cells. Genomic evidence also supports the association of severe COVID-19 with various inflammatory syndromes, such as neuronal COVID-19, acute respiratory distress syndrome, vascular inflammation, and multisystem inflammatory syndrome. In conclusion, this study suggests that gene expression profiling plays a significant role in elucidating the etiology of severe COVID-19.

Transcriptomic analysis and fusion gene identifications of midbody remnants released from colorectal cancer cells reveals they are molecularly distinct from exosomes and microparticles.

Previously, we reported that human primary (SW480) and metastatic (SW620) colorectal (CRC) cells release three classes of membrane-encapsulated extracellular vesicles (EVs); midbody remnants (MBRs), exosomes (Exos), and microparticles (MPs). We reported that MBRs were molecularly distinct at the protein level. To gain further biochemical insights into MBRs, Exos, and MPs and their emerging role in CRC, we performed, and report here, for the first time, a comprehensive transcriptome and long noncoding RNA sequencing analysis and fusion gene identification of these three EV classes using the next-generation RNA sequencing technique. Differential transcript expression analysis revealed that MBRs have a distinct transcriptomic profile compared to Exos and MPs with a high enrichment of mitochondrial transcripts lncRNA/pseudogene transcripts that are predicted to bind to ribonucleoprotein complexes, spliceosome, and RNA/stress granule proteins. A salient finding from this study is a high enrichment of several fusion genes in MBRs compared to Exos, MPs, and cell lysates from their parental cells such as MSH2 (gene encoded DNA mismatch repair protein MSH2). This suggests potential EV-liquid biopsy targets for cancer detection. Importantly, the expression of cancer progression-related transcripts found in EV classes derived from SW480 (EGFR) and SW620 (MET and MACCA1) cell lines reflects their parental cell types. Our study is the report of RNA and fusion gene compositions within MBRs (including Exos and MPs) that could have an impact on EV functionality in cancer progression and detection using EV-based RNA/ fusion gene candidates for cancer biomarkers.

Multi-level profiling of the Fmr1 KO rat unveils altered behavioral traits along with aberrant glutamatergic function

Fragile X syndrome (FXS) is the most common cause of inherited intellectual disabilities and the most prevalent monogenic cause of autism. Although the knockout (KO) of the Fmr1 gene homolog in mice is primarily used for elucidating the neurobiological substrate of FXS, there is limited association of the experimental data with the pathophysiological condition in humans. The use of Fmr1 KO rats offers additional translational validity in this regard. Therefore, we employed a multi-level approach to study the behavioral profile and the glutamatergic and GABAergic neurotransmission status in pathophysiology-associated brain structures of Fmr1 KO rats, including the recordings of evoked and spontaneous field potentials from hippocampal slices, paralleled with next-generation RNA sequencing (RNA-seq). We found that these rats exhibit hyperactivity and cognitive deficits, along with characteristic bidirectional glutamatergic and GABAergic alterations in the prefrontal cortex and the hippocampus. These results are coupled to affected excitability and local inhibitory processes in the hippocampus, along with a specific transcriptional profile, highlighting dysregulated hippocampal network activity in KO rats. Overall, our data provide novel insights concerning the biobehavioral profile of FmR1 KO rats and translationally upscales our understanding on pathophysiology and symptomatology of FXS syndrome.

Pan-cancer Genomic Analysis of AXL Mutations Reveals a Novel, Recurrent, Functionally Activating AXL W451C Alteration Specific to Myxofibrosarcoma.

Myxofibrosarcoma (MFS) is a common soft tissue sarcoma of the elderly that typically shows low tumor mutational burden, with mutations in TP53 and in genes associated with cell cycle checkpoints (RB1, CDKN2A). Unfortunately, no alterations or markers specific to MFS have been identified and, as a consequence, there are no effective targeted therapies. The receptor tyrosine kinase AXL, which drives cellular proliferation, is targetable by new antibody-based therapeutics. Expression of AXL messenger RNA is elevated in a variety of sarcoma types, with the highest levels reported in MFS, but the pathogenic significance of this finding remains unknown. To assess a role for AXL abnormalities in MFS, we undertook a search for AXL genomic alterations in a comprehensive genomic profiling database of 463,546 unique tumors (including 19,879 sarcomas, of which 315 were MFS) interrogated by targeted next-generation DNA and/or RNA sequencing. Notably, the only genomic alterations recurrent in a specific sarcoma subtype were AXL W451C (n = 8) and AXL W450C (n = 2) mutations. The tumors involved predominantly older adults (age: 44 to 81 [median: 72] y) and histologically showed epithelioid and spindle-shaped cells in a variably myxoid stroma, with 6 cases diagnosed as MFS, 3 as undifferentiated pleomorphic sarcoma (UPS), and 1 as low-grade sarcoma. The AXL W451C mutation was not identified in any non-sarcoma malignancy. A review of publicly available data sets revealed a single AXL W451C-mutant case of UPS that clustered with MFS/UPS by methylation profiling. Functional studies revealed a novel activation mechanism: the W451C mutation causes abnormal unregulated dimerization of the AXL receptor tyrosine kinase through disulfide bond formation between pairs of mutant proteins expressing ectopic cysteine residues. This dimerization triggers AXL autophosphorylation and activation of downstream ERK signaling. We further report sarcomas of diverse histologic subtypes with AXL gene amplifications, with the highest frequency of amplification identified in MFS cases without the W451C mutation. In summary, the activating AXL W451C mutation appears highly specific to MFS, with a novel mechanism to drive unregulated signaling. Moreover, AXL gene amplifications and messenger RNA overexpression are far more frequent in MFS than in other sarcoma subtypes. We conclude that these aberrations in AXL are distinct features of MFS and may aid diagnosis, as well as the selection of available targeted therapies.

Biomarker selection depends on gene function and organ: the case of the cytochrome P450 family genes in freshwater fish exposed to chronic pollution.

Pollution and its effects have been of major concern in recent decades. Many strategies and markers have been developed to assess their effects on biota. Cytochrome P450 (CYP) genes have received significant attention in this context because of their relationship with detoxification and activation of exogenous compounds. While their expression has been identified as a pollution exposure biomarker, in most cases, it has been tested only after acute exposures and for CYP genes associated with exogenous compounds. To elucidate CYP gene expression patterns under chronic pollution exposure, we have used the silverside Basilichthys microlepidotus as a model, which inhabits the Maipo River Basin, a freshwater system with different pollution levels. We performed next-generation RNA sequencing of liver and gill tissues from polluted and non-polluted populations. We found most CYP genes were not dysregulated by pollution, and the seven genes that were present and differentially expressed in liver and gill were mainly downregulated. Three CYP genes associated with exogenous compounds showed differential expression in the gill, while four CYP genes associated with endogenous compounds showed differential expression in the liver. The findings presented here highlight the importance of CYP genes, his family, tissues and his interaction in the context of pollution biomarkers use.

Divergent transcriptomic signatures from putative mesenchymal stimuli in glioblastoma cells

AbstractIn glioblastoma, a mesenchymal phenotype is associated with especially poor patient outcomes. Various glioblastoma microenvironmental factors and therapeutic interventions are purported drivers of the mesenchymal transition, but the degree to which these cues promote the same mesenchymal transitions and the uniformity of those transitions, as defined by molecular subtyping systems, is unknown. Here, we investigate this question by analyzing publicly available patient data, surveying commonly measured transcripts for mesenchymal transitions in glioma-initiating cells (GIC), and performing next-generation RNA sequencing of GICs. Analysis of patient tumor data reveals that TGFβ, TNFα, and hypoxia signaling correlate with the mesenchymal subtype more than the proneural subtype. In cultured GICs, the microenvironment-relevant growth factors TGFβ and TNFα and the chemotherapeutic temozolomide promote expression of commonly measured mesenchymal transcripts. However, next-generation RNA sequencing reveals that growth factors and temozolomide broadly promote expression of both mesenchymal and proneural transcripts, in some cases with equal frequency. These results suggest that glioblastoma mesenchymal transitions do not occur as distinctly as in epithelial-derived cancers, at least as determined using common subtyping ontologies and measuring response to growth factors or chemotherapeutics. Further understanding of these issues may identify improved methods for pharmacologically targeting the mesenchymal phenotype in glioblastoma.

Altered N6-Methyladenosine Modification Patterns and Transcript Profiles Contributes to Cognitive Dysfunction in High-Fat Induced Diabetic Mice.

N6-methyladenosine (m6A) constitutes the paramount post-transcriptional modification within eukaryotic mRNA. This modification is subjected to stimulus-dependent regulation within the central nervous system of mammals, being influenced by sensory experiences, learning processes, and injuries. The patterns of m6A methylation within the hippocampal region of diabetes cognitive impairment (DCI) has not been investigated. A DCI model was established by feeding a high-fat diet to C57BL/6J mice. m6A and RNA sequencing was conducted to profile the m6A-tagged transcripts in the hippocampus. Methylated RNA immunoprecipitation with next-generation sequencing and RNA sequencing analyses yielded differentially m6A-modified and expressed genes in the hippocampus of DCI mice, which were enriched in pathways involving synaptic transmission and axonal guidance. Mechanistic analyses revealed a remarkable change in m6A modification levels through alteration of the mRNA expression of m6A methyltransferases (METTL3 and METTL14) and demethylase (FTO) in the hippocampus of DCI mice. We identified a co-mediated specific RNA regulatory strategy that broadens the epigenetic regulatory mechanism of RNA-induced neurodegenerative disorders associated with metabolic and endocrine diseases.

B cell immune response to human bocaviruses.

Human bocaviruses (HBoVs) have been demonstrated in respiratory and gastrointestinal infections; however, the immune response to them has not been studied in detail. In this study, we investigated the B cell immune responses to HBoV1 and HBoV2, representing two different species of bocaviruses in humans. We analyzed the effects of stimulations with HBoV1 and 2 virus-like particles (VLPs) and of co-stimulation with HBoV1-rhinovirus (RV) on cells of the immune system by flow cytometry, transcriptomics, and luminometric immune assays. Human B cells, and particularly B regulatory cells (Breg cells), showed an increased immune response to HBoV1-VLPs stimulation. These immune responses were also supported by increased IL-1RA and PDL1 expressions in IL-10+ B cells from peripheral blood mononuclear cells (PBMCs) stimulated with HBoV1-VLPs. In addition, increased levels of IL-10 and IL-1RA were determined in the supernatants of PBMCs following HBoV1-VLPs stimulation. HBoV1-VLPs and RV co-stimulation increased the IL-10+ B cell population. Transcriptome analysis by next-generation RNA sequencing showed an increased expression of IL-10 signalling and Breg cell markers in PBMCs stimulated with HBoV1-VLPs. Furthermore, TGF-β and chemoattractants MIP-1α, MIP-1β and IP10 protein levels were high in the supernatants of PBMCs stimulated with HBoV1-VLPs. The findings demonstrate that in Breg cells, IL-10 signalling pathways, and anti-inflammatory activity are induced by HBoV1, which can explain the often mild nature of the disease. In addition, the immune regulatory response induced by HBoV1-VLPs may indicate a potential immunomodulatory role of HBoV1 on the immune system and may represent an immune regulatory strategy.

Identifying genetic susceptibility to Aspergillus fumigatus infection using collaborative cross mice and RNA-Seq approach.

Aspergillus fumigatus (Af) is one of the most ubiquitous fungi and its infection potency is suggested to be strongly controlled by the host genetic background. The aim of this study was to search for candidate genes associated with host susceptibility to Aspergillus fumigatus (Af) using an RNAseq approach in CC lines and hepatic gene expression. We studied 31 male mice from 25 CC lines at 8 weeks old; the mice were infected with Af. Liver tissues were extracted from these mice 5 days post-infection, and next-generation RNA-sequencing (RNAseq) was performed. The GENE-E analysis platform was used to generate a clustered heat map matrix. Significant variation in body weight changes between CC lines was observed. Hepatic gene expression revealed 12 top prioritized candidate genes differentially expressed in resistant versus susceptible mice based on body weight changes. Interestingly, three candidate genes are located within genomic intervals of the previously mapped quantitative trait loci (QTL), including Gm16270 and Stox1 on chromosome 10 and Gm11033 on chromosome 8. Our findings emphasize the CC mouse model's power in fine mapping the genetic components underlying susceptibility towards Af. As a next step, eQTL analysis will be performed for our RNA-Seq data. Suggested candidate genes from our study will be further assessed with a human cohort with aspergillosis.

Canonical Wnt signaling pathway (WSP) alterations in metastatic prostate cancer.

221 Background: WSP aberrations have been implicated in the progression and treatment resistance of metastatic prostate cancer (mPCa). Sustained exposure of PCa to androgen deprivation has been shown to regulate the canonical WSP, induce beta-catenin expression and stabilization, and promote androgen independence. This study interrogates a large dataset to characterize the molecular and immune landscape of canonical WSP-altered mPCa and associations with survival outcomes. Methods: We examined PCa patients who underwent tissue-based DNA and RNA sequencing through a CLIA-certified commercial assay (Caris Life Sciences). Comprehensive genomic profiling included next-generation sequencing (NGS) of DNA (592-gene panel or whole exome) and RNA (whole transcriptome), on tissue samples from primary prostatic and/or metastatic sites. Tumor samples were categorized as having aberrant canonical Wnt signaling (WSP-act) if they harbored activating mutations in CTNNB, pathogenic RSPO2 fusions, or inactivating mutations in either APC or RNF43. Subset analysis was conducted in mPCa samples with microsatellite stability (MSS), after excluding those with RNF43 (p.G659fs*) mutations which are passenger events in MSI-high cancers. Statistical significance was determined using Fisher's Exact or X2 tests with Benjamini–Hochberg correction. Overall survival (OS) data was obtained from insurance claims and was analyzed using Kaplan-Meier estimation. Results: Of all PCa samples (4,150) (N = 4,150), 17.4% (n = 722) were classified as WSP-act; 299 (41.4%) were from the primary tumor and 423 (58.6%) were from metastatic sites. WSP-act tumors were more prevalent in brain (46%), liver (31%), lung (31%) and bone (25%) metastases compared to primary prostate cancers (10%). The MSS mPCa subgroup (N=1,638) comprised 321 WSP-act and 1,317 WSP-WT tumors. WSP-act mPCa exhibited more SPOP mutations (14.2% vs. 6.3%, p < 0.001) as well as higher mRNA expression of CST1, NKD1, and ODAM (known canonical WSP activators) compared to WSP-WT cases. M2 macrophages predominated the tumor microenvironment, with a higher median fraction in WSP-act tumors compared to WSP-WT (7.2 vs. 6.1%, p < 0.001). ROR1 gene expression was significantly elevated in WSP-act tumors from both primary and metastatic sites. There was no significant difference in OS outcomes between WSP-act and WSP-WT mPCa patients. Conclusions: WSP-act prostate cancer demonstrated a pronounced upregulation of ROR1 gene expression, underscoring its potential involvement in the crosstalk between the canonical and non-canonical Wnt signaling pathways. Additionally, the augmented levels of M2 macrophages in WSP-act tumors, combined with the reported role of ROR1 in tumor immunosuppression, suggests that ROR1 may contribute to immune evasion in WSP-act mPCa. Our findings may provide rationale for developing novel therapeutic strategies targeting Wnt-activated prostate cancers.